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Abstract 

This work presents analytical formulas derived for evaluating vapor transport of a volatile solvent for 
an isolated multicomponent droplet in a quiescent environment, based on quasi-steady-state approximation. 
Among multiple solvent components, only one component is considered to be much more volatile than the 
rest such that other components are assumed to be nonvolatile remaining unchanged in the droplet during 
the process of (single-component) volatile solvent evaporation or condensation. For evaporating droplet, 
the droplet size often initially decreases following the familiar “cP law” at an accelerated rate. But toward 
the end, the rate of droplet size change diminishes due to the presence of nonvolatile cosolvent. Such 
an acceleration-deceleration reversal behavior is unique for evaporating multicomponent droplet, while the 
droplet of pure solvent has an accelerated rate of size change all the way through the end. This reversal 
behavior is also reflected in the droplet surface temperature evolution as “S-shaped” curves. However, 
a closer mathematical examination of conditions for acceleration-deceleration reversal indicates that the 
acceleration phase may disappear when the amount of nonvolatile cosolvent is relatively small and ambient 
vapor pressure is relatively high. Because the net effect of adding nonvolatile cosolvent is to reduce the 
mole fraction of the volatile solvent such that the saturation vapor pressure is lowered, vapor condensation 
onto the multicomponent droplet is predicted to occur when the ambient vapor pressure is subsaturated with 
respect to that for the pure volatile solvent. In this case, the droplet will grow asymptotically toward a finite 
size. But when the ambient vapor pressure becomes supersaturated with respect to that for the pure volatile 
solvent, the condensation growth of droplet can continue indefinitely without bound. 


1 Introduction 

In Aerosol Jet® additive manufacturing, the ink, consisting of both volatile and nonvolatile solvents besides 
the solids, is atomized into a dense mist of fine aerosol droplets usually about 50 nL/cc or more. The ink 
mist, with droplets of typical diameters in the range between 1 to 5 /xm, is then carried through conditioning 
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devices to an aerodynamic focusing nozzle for precision material deposition onto a substrate (cf. Zollmer 


et al. 2006| Kahn| 2007} Christenson et al. 20111. The primary functionality of the volatile component i 


IS 


to reduce ink viscosity so as to enable effective atomization. Once atomized, the ink droplets are desired 
to have high solid content especially when they arrive at the substrate. Higher solid content corresponds to 
higher viscosity as preferred for retaining the printed features at high mass deposition rate, which may not 
be easily achieved with an ink too “runny” on the substrate. Therefore, removal of the volatile component in 
ink droplets before deposition plays an important role in high-quality Aerosol Jet® printing. In fact, the ink 
for Aerosol Jet® printing is typically formulated to contain solvents with one volatile component and one 
non-volatile component. The reason for having a non-volatile component is, while the volatile component 
is mostly removed, to retain cohesivity of the printed features for further solidihcation processing such as 
controlled drying, sintering, with appropriate annealing. 

Because the ink droplets in the Aerosol Jet® printing process travel through a transportation channel 
from the atomizer to the deposition nozzle before arriving the substrate, opportunities are available for “in¬ 
flight” droplet solid content adjustment via the channel shape design, selection of mist flow rate and sheath 
flow rate, variation of the ratio between volatile and nonvolatile solvent components in the ink formulation, 
etc. Thus, non-contact deposition of ink droplets with high solid content becomes possible with the Aerosol 
Jet® technology, enabling single-pass high-aspect-ratio feature printing. When designing-optimizing an 
Aerosol Jet® printer, formulating an ink, or developing a process recipe for a particular application, under¬ 
standing the physical principles behind volatile solvent transport for a multicomponent ink droplet becomes 
practically important with an analytical theory that can provide reasonable estimate of relevant time scales 
involved therewith. 

The main purpose of the present work is to derive analytical formulas for conveniently evaluating the 
evaporation as well as condensation characteristics of an isolated ink droplet. It can be easily used to provide 
theoretical guidance to ink formulation as well as Aerosol Jet® printing recipe development. For consis¬ 
tent printing performance, relative motions of ink droplets in the Aerosol Jet® mist stream are minimized 
by design as the mist flowing through the deposition nozzle and exposing to a coflowing sheath gas with 
controllable solvent vapor contents, where substantial droplet evaporation occurs. Thus, a model based on 
vapor transport of a volatile solvent for an isolated multicomponent droplet in a quiescent surrounding gas 
is expected to be quite relevant. 

In a bigger picture, the evaporation behavior of aerosol droplets is of interested to a variety of applica- 


tions such as spray drying (|L 
aerosol drug delivery to the 1 
name a few. It is a subject stu 

n and Gentry| 20031, liquid fuel combustion (|Aggarwal[ 19981, pharmaceutical 
ungs (Peng et al.| 2000|l, atmospheric science ([Pruppacher and Klettl 19781, to 

died theoretically by many authors (e.g.. 

Maxwell! 1^78[ Pruppacher and Klett 

19781 Ravindran and Davis| 

1982[ Kulmala et al.| 1993[ Widmann and Davis 1997||. Analytical formulas 

were presented for evaporation of a single component (unitary) droplet (Maxwell 1878[ Pruppacher and 


Klett[ 19781 or evaporation of a two-component (binary) droplet with constant and uniform composition 


(Kulmala et al. 1993| Vesala et al. 19971. In mathematical derivations, the assumptions of “quasi-steady- 
state” and “diffusion-controlled” mass transfer have typically been used. Coupled nonlinear equations for 
multicomponent evaporation of a multicomponent droplet were derived with solutions computed numeri¬ 
cally (Ravindran and Davis 1982| Widmann and PavTsl 1997| Wilms 20051. Comparison between theory 
and experiment have been presented by Wilms ( |2005[ ); Davies et al. (20121, among others. 

For Aerosol Jet® printing, the in-process evaporation of only a single (volatile) solvent component from 
a multicomponent droplet is of interest for most of the times. Although the involvement of multiple sol¬ 
vent components in the ink droplets makes mathematical derivations somewhat complicated, an analytical 
formula can still be obtained and can be conveniently used for effective technology development. In what 
follows, detailed mathematical derivation steps with simplifying assumptions are presented in section 2, 
followed by results and discussion in section 3 for both evaporating droplet and growing droplet by conden¬ 
sation. Finally, the key predictions are summarized in section 4. 
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2 Mathematical Derivation 


A spherical multicomponent liquid droplet is considered here in a quiescent surrounding gas, with one of 
its solvent components (the volatile component) evaporating (or condensing) while the amount of other 
(nonvolatile) components remaining unchanged during the time scale of interest. The ink droplets suitable 
for Aerosol Jet® printing typically have diameters in the range of 1 to 5 /rm, containing about 10% volume 
of solids (whose only role is considered to passively occupy some of the droplet volume). Thus, the size 
of an ink droplet after complete removal of the volatile solvent should be about half of its original size, 
still much larger than the mean free path of the surrounding gas molecules under ambient condition. It is 


therefore not unreasonable to ignore the Knudsen number effect. Similar to Feng (20131, the spherically 
symmetric mass transfer of the volatile solvent vapor in an incompressible continuous gas medium with 


density pg and diffusivity D is governed by (Bird et ah) 1960 p. 557) 


dpv 


+ V 


dpv 

dr 


= ^^(d^] {R<r<oo) 


dr 


( 1 ) 


where p^ denotes the mass concentration of the solvent vapor, t the time, and r the radial distance from 
the center of the droplet of radius R{i). To satisfy the equation of continuity while accounting for the 
effect of volume change during the solvent phase change, the radial velocity v of spherically symmetric 
incompressible flow is simply (Scriven 1959|) 


?; = 1 - 



( 2 ) 


where pi denotes the liquid density of the solvent under consideration. 

At the droplet surface, the mass balance based on Pick’s first law of binary diffusion accounting for the 
bulk flow effect with the gas medium flux being ignored (Bird et ah] 1960[ p. 502) leads to 


dR _ PgD 

di PliPg-Pv,R) 


dpv 

dr 


(3) 


r=R 


where denotes the vapor density at droplet surface (r = R). 

For most volatile solvents used in Aerosol Jet® inks, the density of solvent vapor is much less than 
that of gas medium, i.e., py^R « Pg the diffusivity D may be treated as a constant. Thus, in (|Tjl the 
convection term 



PgD 


PliPg - Pv,r) 


dpv 

dr 


r=R 


dpv 

dr 


Pv,Rj^Rdp> 

Pg 


dr 


may be assumed negligible comparing with other terms. Then, with the terms associated to Pv,r/P g ignored, 
the solution to ([T]l leads to mass balance at the droplet surface 


dR _ Di^pyR Pv,oo) 1 

di 


Pi 



(4) 


wherep^oo denotes the ambient vapor density at r —)• oo. Although exact solution to ( [^ ma y be derived 
(cf., Feng 20131, retaining the complexity due to the transient diffusion flux term R/wnDt seems to be 
unnecessary because significant change of R occurs on a time scale ~ piR^/{-KDpy^R) (with Rq denoting 
the droplet initial radius used here as the characteristic length) which is much greater than K^KttD) for 
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Pi » Pv,R liquid density is much greater than gas density). Hence, the quasi-stationary mass balance 
equation Q may be reasonably approximated by a quasi-steady-state equation 

dR'^ _ 2 D{p^^R - 5O0) 

dt Pi 


which has the same form as standard expression based on Maxwell continuum transport approximation (e.g., 


vapor transport to and from a droplet also involves 
substantial heat transfer owing to the release or absorption of heat of phase change, namely the latent heat. 
The resulting temperature difference between the droplet and its surrounding gas medium causes a flow of 
sensible heat by thermal diffusion. With the quasi-steady-state approximation, energy balance at the droplet 
surface leads to 


Pruppacher and Klett| 1978[ Widmann and DavTs[ 1997 


As expected (Pruppacher and Klett 1978[ p. 417) 


iPv,R Pv,oo) — '^Ft) 


( 6 ) 


where k is the thermal conductivity of the surrounding gas medium, L the latent heat due to volatile solvent 
phase change. Too und Tr the temperature values at r —00 and r = R. Moreover, it is usually reasonable to 
take p^^R as the local saturated vapor density at the droplet surface Psv,R (with “sv” in the subscript denoting 
saturated vapor). Making use of the Clausius-Clapeyron relation, for characterizing a discontinuous phase 
transi tion between liquid and gas, yields an approximate equation (as justified in [Pruppacher and Klett[ 
p. 420) relating local saturated vapor density difference between r = R and r —)> 00 to {Tr — Too) as 
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(.Psv,R Psv,oo') — 


PsVyi 




L 


2 rn 
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where denotes the (specific) gas constant for the volatile solvent vapor (which can be calculated as 
the universal gas constant, 8.3143 x 10^ J deg“^ kmol“^, divided by the solvent molecular weight in kg 
kmol“^.) 

Combination of (|^ and ([7]l leads to 


{Psv,R Pv,oo) 



Pv,c 


Psv,c 


1 LD 

psv^oo kToQ 




( 8 ) 


Thus, expressing (p„,_R — Pv,oo) in terms of corresponding vapor pressure Pv,oo = Pv,oo/{^vToo), Psv,oo = 
Psv,oo/i^vToo), and Too yields a more convenient form of ([^ as 


dR‘^ 

dt 


2 Pv,oo \ 

V Psv,oo J 


Pl^yToo piLD 

Psv,00 kT^ 




( 9 ) 


where all the variables become dimensionless with R{t) measured in units of the droplet initial (dimen¬ 
sional) radius Rq and t in units of Rq/{'kD), i.e., R = R/Rq and t = It^D/Rq. Clearly, R{t) decreases 
with t when Pv,oo/Psv,oo < 1 as the solvent is evaporating from the droplet. If the droplet is in a supersat¬ 
urated environment when Pv,oo/Psv,oo > l^ R{t) would increase with t because in this case the solvent is 
condensing onto the droplet. 

If the droplet of multicomponent solution contains only one volatile component, as Aerosol Jet® inks 
are typically formulated, the saturated vapor pressure Psv,oo for tho volatile solvent becomes a function of 
its mole fraction x in the solution. Assuming the multicomponent mixture is an ideal solution such that 
Raoult’s law is applicable, we have 

Psv,00 - xp*sv,00 1 ( 10 ) 

with p*sv,oo denoting the saturated vapor pressure of the pure volatile solvent. Substituting ( jlOj ) to 0 incurs 
complexity due to the fact that x is a function of R{t). With the (dimensionless) volume of nonvolatile 
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component of liquid solution in the ink droplet denoted by 47ri2^/3 and the total volume remaining in 
the droplet when the volatile component is completely removed—the remnant volume—by AttR^ /3 (which 
includes both solids and nonvolatile component of liquid solution), the volume of volatile solvent component 
should be 47r(i?^ — i?^)/3 and 


1 PnMy 

x{t) PvMn R{t)^ - R^ 


( 11 ) 


where pn and p^, and My denote the liquid solvent densities, molecular weights of the nonvolatile and 
volatile components, respectively. Clearly, the mole fraction x is generally a function of time t during the 
volatile solvent evaporation, because R{t) changes with time while Rn and Ry remaining as constants. 
Incorporating ([T0|) and ([TT]) in Q yields 
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2 1 Pv,oo/P^s 
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An analytical solution to ([T^, though apparently complicated, can be written as 


( 12 ) 


(13) 

(14) 

(15) 


j^2 _ 2t(l Pv,oo/P*sv,oo) ^ 
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3AR, 


ttA 
In 


{l-RefjR^-Rl) 
l{R-Re)Hl-Rl) \ 


(16) 


Unfortunately, ( [T^ is in an implicit form of i? as a function t. However, t can be expressed as an explicit 
function of R, which make it straightforward to calculate the value of t from a given value of R for generating 
curves of R versus t. 


3 Results and Discussion 


Without going to the solution, an examination of ( [I^ suggests the possible situation for dR?/dt —)■ 0 as 
R ^ Re with Re denoting the equilibrium droplet radius as defined in ( [T3] ) for the droplet at thermody¬ 
namic equilibrium with ambient vapor pressure py^oo- Of course, this can only happen in the presence of 
the nonvolatile component, i.e., when Rn is nonzero such that B is nonzero. In the absence of the non¬ 
volatile cosolvents (i.e., Rn = 0) that can influence the equilibrium vapor pressure of the volatile solvent, 
B disappears and R^ simply becomes a linear function of t for a droplet containing a single solvent. 

While (nondimensional) Rn and Ry are subjected to constraint as Rn < Rr < I (measured in units of 
droplet initial radius Rq), there is no upper limit to the equilibrium radius Re. According to ( [T3] ), however, 
i?e > 1 can only happen when R/^py^^o / 0 in a subsaturated environment, i.e., Py^oo < P*sv,oo, while 
satisfying the condition p /Pv,oo < l/x(0) where 


^ I PnAdyRn 

x(0) PvMn{l - R/:) 


(17) 
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Table 1: Solvent parameter values under standard ambient condition for density pi (kg m“^), molecular 
weight Mw (kg kmol“^), saturated vapor pressure p*sv,oc (Pa), latent heat of vaporization L (J kg“^), 
diffusivity in air D (m^ s“^), and dimensionless A defined in ( [l4| ) using the thermal conductivity in air 
k = 0.024 (J m-i s"^ K"!). 


Solvent 

Pi 

My, 

P*sv,oo 

L 

D 

A 

water 

1000 

18.02 

2.99 X 10^ 

2.26 X 10® 

2.50 X 10"® 

1.67 X 10® 

butanol 

810 

74.12 

5.53 X 102 

7.08 X 10® 

8.00 X 10“® 

6.21 X 10^ 

EG 

1113 

62.07 

8.00 

— 

— 

— 


In the case of Pv,oo = a;(0)p*s.i;^oo^ we have iie = 1 and hence dR^/dt = 0 forever because the ambient 
vapor pressure equals the saturation value starting from t = 0. Volatile solvent evaporation occurs when 
Re < 1, which makes R{t) (> Re) decreases with time until R{t) —)• Re. Vapor condensation onto the 
droplet can occur either when Re > 1 such that R{t) (< Re) increases with time until reaching Re, or when 
Pv,cx) > P*sv,oo- In a supersaturated environment, namely, pv,cx) > P*sv,cx), we can have Re < 0 < Rr < I 
according to ( jlsj ), in which case the negative Re loses its physical meaning, but dR?/dt > 0 as the droplet 
is growing due to Pv,oc/p*sv,oc > 1 as indicated by ( [T^ . Supersaturated vapor will condense onto any 
surfaces, such as the channel walls, besides the surface of a multicomponent droplet. Noteworthy here is 
that ( [T3] ) also suggests an indefinte growth of droplet (without bound) as Pv,oo P*sv,oo^ because it leads 
to Re oo. 

For illustrations in this section, water and butanol are used as realistic examples for the volatile solvent, 
and ethylene glycol (EG) for the nonvolatile solvent. They have actually been used in several Aerosol 
Jet® ink formulations. Table 1 shows the typical parameters for those solvents, where only the liquid 
solvent density pi and molecular weight are relevant for the nonvolatile ethylene glycol. The ambient 
temperature Too is assumed to be 300 K and ambient pressure to be 10^ Pa. Thus, the values of hjiM^Tea) 
for water and butanol are 16.3 and 21.0, respectively. As a consequence, the value of Psv,R is generally 
less than psv,oo according to ([^. The values of Pv,r/ Pg in (j^, which is less than Psv,oo/^^^, are therefore 
< 0.03 and 0.0055 for water and butanol, as a posteriori verification of negligible terms associated with 
factor Pv,r/ Pg in our mathematical derivation. 

From (16 1 it is clear that the effects of nonvolatile cosolvent are proportional to the value of B j A. While 
the value of A remains constant for a given solvent, the value of B, according to ( [T5] ), is proportional to 
the volume fraction of nonvolatile cosolvent R^ (among other factors such as the volatile solvent molecular 
weight My, etc.) Thus, a simple inspection of the physical property values in Table 1 would suggest that 
a butanol-EG droplet should exhibit much more significant effects of nonvolatile solvents than a water-EG 
droplet. 


3.1 Evaporating Droplet 

When the equilibrium radius Re < 1, the volatile solvent will evaporate from the droplet of initial R = 1 
until R —)■ Re- In the absence of the nonvolatile cosolvent, i.e., = 0, we have B = 0 and ( [T^ recovers 

the familiar classical “d^ law” for evaporation of a single component droplet with a constant dR?/dt. In this 
case, the volatile solvent can be completely removed in a finite time 


to 


■kA{ 1 — Ry) 
2(1 Pv,oo /P*sv 


(18) 
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Figure 1: Evaporating droplet radius evolution with time for water-EG and butanol-EG mixtures, with 
ethylene glyeol (EG) serving as the nonvolatile cosolvent with a volume fraction of = 0.4. The remnant 
volume of droplet (excluding that of the volatile solvent) is denoted as AirR^/S here in this figure with the 
remnant volume fraction R^ = 0.5. The solid curve and dashed curve correspond to Pv,oo/p*sv,oo = 0 and 
0.2, with = 1.23 and 1.71 for water-EG droplet while = 4.81 and 5.49 for butanol-EG droplet. 
The dotted curve is for pure volatile solvent with Rf^ = 0 and Rf. = 0.5 at Pv,oo/p*sv,oo = 0, as a reference. 


A nonzero Rn leads to a nonzero B that causes a fundamental change in the evaporation behavior: the 
volatile solvent cannot be completely removed in a finite time because now dR?/dt —)• 0 as ii —)• i?e due to 
diminishing mole fraction of the volatile solvent x in the droplet. As a consequence, it is expected to take 
forever to remove the last bit of the volatile solvent from a droplet in the presence of a nonvolatile cosolvent. 

If the droplet has a volume fraction of nonvolatile cosolvent R^ = 0.4 and remnant volume fraction 
Rf = 0.5 (namely, Rr = 0.794 as the remnant radius), curves of R{t) are shown in Eig. 1 for evaporating 
droplets of water-EG mixture and butanol-EG mixture with 10% solids and 40% nonvolatile solvents in 
the initial volume, at Pv,oo/p*sv,oo = 0 (solid curve) and 0.2 (dashed curve) which is about 25% of x(0) 
(« 0.795) for the water-EG droplet and nearly 50% of x(0) (ss 0.432) for the butanol-EG droplet. As a 
reference, also shown in Eig. 1 is the dotted curve representing the case for Rn = 0 and Rf = 0.5 at 
Pv,oo/p*sv,oo = 0 (according to the “d^ law”) for a single solvent droplet. Eor convenience of displaying 
results, t* in all figures here is used to denote a normalized time with the value of to in (181 at Pv,oc /p*sv,oc = 
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0 such as 


t* = 


2t 


ttA{1-R^) ■ 

Accordingly, the dimensional time corresponding to t* = 1 can be calculated as 

A{1 - Rl)Rl 
ZD ’ 


(19) 


( 20 ) 


with i?o and D taking appropriate dimensional values. For example, a water droplet of Rq = 10“® m (one 
micron radius, as the typical mean ink droplet radius in the Aerosol Jet® mist) containing 50% volume 
fraction of solids (i.e., a remnant volume fraction of R^ = 0.5) will completely evaporate (at t* = 1) in 1.24 


milliseconds, according to (201 with parameters in Table 1. But for a butanol droplet of the same size with 
the same volume fraction of solids, the time for complete evaporation will be 1.43 milliseconds. Increasing 
the droplet size by a factor of 2 will increase the dimensional time by a factor of 4 for droplets of the same 
composition, as indicated by (201. If the remnant volume fraction R^ is reduced to 0.2 (or 0.1), the time for 


complete evaporation of pure solvent-solid droplets would increase by a factor of 1.78 (or 2.12) from that 
for R^ = 0.5. 

Strictly speaking, removal of all the volatile solvent from a multicomponent droplet containing non¬ 
volatile cosolvent would take forever, i.e., t* —)• oo. However, it is still of practical interest to evaluate the 
time for “almost” complete evaporation. Therefore, a time denoted as for R{t) to reach 1.01 x R^ is 
used here for evaluating the (finite) time t* needed for reaching R k, R^,. The value of can be calculated 
by substituting 1.01 x Rg for R in (161. For example, the values of for a water-EG droplet with non¬ 
volatile cosolvent volume fraction R^ — 0-4 and remnant volume fraction R^ = 0.5 at Pv,oo/p*sv,oo = 0 
and 0.2 are 1.23 and 1.71, respectively. But for a butanol-EG droplet with 72^ = 0.4 and Rf = 0.5 at 
Pv,oo/p*sv,oo = 0 and 0.2, the values of become are 4.81 and 5.49, respectively. Thus, it takes a water- 
EG droplet in Eig. 1 with an initial radius of one micron (Rq = 10“® m) about (t*Q]^ x 1.24) 1.53 ms and 
2.12 ms to “almost” complete the evaporation process at Pv,oo/p*sv,oo = 0 and 0.2, whereas a butanol-EG 
droplet (of Rq = 10“® m) in Eig. 1 needs about x 1.43) 6.88 ms and 7.85 ms. Obviously, the presence 
of the nonvolatile cosolvent (with a volume fraction of Rf^ = 0.4) hinders the volatile solvent evaporation 
of a butanol-EG droplet much more significantly than that of a water-EG droplet. This is expected from the 
fact that the value of B jA for the butanol-EG droplet (= 0.5217) is much greater than that for the water-EG 
droplet (= 0.0358), at Pv,oolp*sv,oo = 0. In the case of Pv,oo/p*sv,oo = 0.2, the values of the equilibrium 
radius Rg become 0.810 for the water-EG droplet and 0.872 for the butanol-EG droplet, respectively. 

If the volume fraction of nonvolatile solvent is reduced to R^ = 0.1 with the remnant volume fraction 
Rf = 0.2, Eig. 2 shows that R{t) for the water-EG droplet almost coincides with that for Rn = 0 at 
Pv,oo/p*sv,oo = 0 in a grant scale, until R{t) becomes very close to Re (= 0.585). Although having 
Pv,oo/p*sv,oo = 0.2 does not cause much change in the value of Re (becoming 0.593 with x(0) = 0.96), 
it noticeably reduces the evaporation rate for the water-EG droplet. Despite the substantial reduction of Rn 
as well as Rr, the butanol-EG droplet still exhibits significant nonvolatile cosolvent effect. Eess nonvolatile 
solvent fraction leads to larger x(0) (= 0.83). Thus, Pv,oo/p*sv,oo = 0.2 corresponds to < 25% of vapor 
equilibrium, which is comparable to the case of water-EG droplet shown in Eig. 1. However, Eig. 2 
indicates more significant nonvolatile cosolvent effect than that for the water-EG droplet in Eig. 1, because 
the value of B/A in this case is 0.130 (still a few times larger than 0.0358 for water-EG in Eig. 1), at 
Pv,oo/p*sv,oo = 0. Eor the case of Pv,oo/p*sv,oo = 0.2, Re of the butanol-EG droplet becomes 0.622 as 
the final equilibrium droplef radius. Wifh R^ = 0.1 and R^ = 0.2 af Pv,oo/p*sv,oo = 0 and 0.2, fhe 
values of fjg^ for a wafer-EG droplef are 1.06 and 1.37, respectively. Buf for a bufanol-EG droplef fhey 
become are 1.94 and 2.57, respectively. Thus, if fakes a wafer-EG droplef in Eig. 2 wifh an inifial radius 
of one micron (Rq = 10“® m) abouf RIqi x 1.24 x 1.78) 2.34 ms and 3.02 ms fo “almosf” complefe fhe 
evaporation process af Pv,oo/p*sv,oo = 0 and 0.2, whereas a bufanol-EG droplef (of Rq = 10“® m) in Eig. 2 








t* 



Figure 2: As in Fig. 1 but for Rf^ = 0.1 and = 0.2, at Pv,cx/p*sv,oo = 0 and 0.2, with = 1.06 and 
1.37 for water-EG droplet while = 1.94 and 2.57 for butanol-EG droplet. 


needs about x 1.43 x 1.78) 4.94 ms and 6.54 ms. Interestingly, the reduction of nonvolatile cosolvent 
volume fraction causes a decrease of time needed for almost complete evaporation of the butanol-EG droplet 
(comparing that in Eig. 1), while the water-EG droplet needs more time, due to intricate effects of the factor 
of 1.78 as a consequence of (1 — with a reduction of the value of for smaller nonvalotile cosolvent 
volume fraction 7?^. 

As an apparent general trend shown in Eig.l and more noticeably so in Eig. 2, the radius of evaporating 
multicomponent droplet tends to follow the law” with increasing magnitude of slope initially when 
7? ~ 1, and then has the magnitude of slope decrease gradually with dR/dt —)• 0 as 77 —)• TZe- In other 
words, evaporating droplet size change is accelerated in the beginning and then decelerated toward the end, 
if a nonvolatile cosolvent is present. Such an acceleration-deceleration reversal behavior is unique in the 
presence of nonvolatile cosolvent in the multicomponent droplet, as also noted by previous authors with 
both experiments and numerical computations (e.g., Widmann and Davis] 1997t Wilms 20051. In contrast, 
the d^-law behavior of accelerating size reduction all the way to 77 = 77e is clearly illustrated by the dotted 
curves for pure solvent cases, which is more obvious in Eig. 2 for the case of R\ = 0.1. 

To have a closer examination of the conditions for acceleration-deceleration reversal phenomenon, it is 
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straightforward to derive the formula for d^R/dr from (|12[) as 


d^R _ (1 - Pv,oo/P*sv,oof{R^ - Rl) 


A {R^ -Rif - 2B {R^ - Rl) - 3BRl 


dt^ 7r2i?2 (^3 _ ^3) + Bf 

which indicates an inflexion point {d^R/dt^ = 0) for R^ > Rl at R = Ri with 


( 21 ) 


O O B 
Ri ~ ^ 


1 1 + ^ _ 1 


B 


( 22 ) 


Thus, we have d?R/dt^ < 0 for i? > i?* and d^R/dt^ > 0 for i? < i?j with d?R/dt^ —)• 0 as i? —)• i?e 
In an explicit form (22) suggests that the value of Ri increases with B which increases with Pv,oo/p*sv,oo 


among others. It is theoretically possible for Ri to become greater than 1. When that happens, the droplet 
evaporation process would not exhibit the acceleration-deceleration reversal behavior; the droplet radius 
reduction process should show a decelerating rate all the way to the end. An evaluation of Ri according to 
(22 1 shows that = 0.887 and 0.929 at = 0 and 0.2 for the water-EG droplet in Fig. 1, which 

indeed exhibits the acceleration-deceleration reversal behavior (because Ri < 1). But for the butanol-EG 
droplet in Fig. 1, which do not seems to show obvious accleration-deceleration reversal, the values of Ri 
are 1.002 and 1.101 (both are greater than 1) at Pv,oo/p*sv,oc = 0 and 0.2. With a reduction of B via 
reduction of the nonvolatile cosolvent volume fraction R^ to 0.1 as in Fig. 2 (for Rf = 0.2, where the 
acceleration-deceleration reversal is observable), the butanol-EG droplet would have Ri = 0.723 and 0.773 

a-tPv,oo/p*sv,oo = 0 and 0.2. 


3.2 Condensational Growth of Droplet 

Droplet can grow by condensation of the volatile solvent vapor when Pv,oo/p*sv,oo > a;(0), which corre¬ 
sponds to i2e > 1 if Pv,oo/p*sv,oo < again until R reaches Rg. According to ( [T3| ), R^ ^ oo at the 
limit of Pv,oo/p*sv,oc = suggesting that the droplet can grow indefinitely without bound. However, for 
Pv,oo/p*sv,oo > the negative value of Re loses physical meaning, as it monotonically increases from 
—oo with pv,oc increasing from p*sv,oo- It should be noted that ( [13] ) can be rewritten as Pv,oo/p*sv,oc 
= {1 -|- pnMyRl/[pyMn{Rl — 7?^)]}“^, which can only relate to the volatile solvent mole fraction x when 
Pv,oo/p*sv,oo < 1, corresponding to Re > Rr- In that case, the droplet size will change until R = Re 
when the volatile solvent mole fraction x{t) in the droplet reaches the thermodynamic equilibrium value of 
Pv,oo/ P*sv,oo- 

The condensation growth behavior of a butanol-EG droplet is illustrated in Fig. 3, with nonvolatile 
cosolvent volume fraction R^ = 0.4 and remnant volume fraction Rl = 0.5, at Pv,oo/p*sv,oo = 0.6, 0.9, 
0.99, 1.01, and 1.10. Because the value of x(0) is 0.432 for Rl = 0.4 and Rl = 0.5, condensation 
growth will occur whenever Pv,oo/p*sv,oo > 0.432. The values of Re corresponding to Pv,oo/p*sv,oo = 
0.6, 0.9, and 0.99 are 1.14, 1.86, and 4.03, respectively. Thus, when Pv,oo/p*sv,oo < condensation 
growth of a multicomponent droplet will practically reach a finite size Re (as t* —)■ cx)). In other words, a 
multicomponent droplet cannot grow indefinitely without bound in an environment of Pv,oo/p*sv,oo < 1 by 
vapor condensation. Useful insights may be gained, however, by examining the normalized time fggg for R 
to reach 0.99 x Rr, similar to for evaporating droplet. For example, with Rl = 0.4 and Rl = 0.5, the 
values of tggg are 13,172, and 8131, for pr^oo/p*sv,oo = 0.6, 0.9, and 0.99, respectively. Thus, a butanol-EG 
droplet of Rq = 10“® m will take 0.029, 0.385, and 18.295 seconds to reach 0.99 x Rr for pr^oo/p*sv,oo = 
0.6, 0.9, and 0.99, given that t* = 1 corresponding to 2.25 milliseconds. When Pv,oo/p*sv,oo > a 
multicomponent droplet will grow indefinitely with time, as seen in Fig. 3 for Pv,oo/p*sv,oo = 1-01 and 
1.10. In this case, the normalized time for 72 = 5 might be a useful measure as the time needed for a 
droplet to grow 125 times of its initial volume. For a butanol-EG droplet with Rl = 0.4 and Rl = 0.5, the 
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Figure 3: As in Fig. 1 (with nonvolatile eosolvent volume fraetion = 0.4 and remnant volume fraetion 
Rf = 0.5) but only for butanol-EG droplet at Pv,oo/p*sv,oo = 0.6, 0.9, 0.99, 1.01, and 1.10 to illustrate 
eharaeteristies of eondensation growth of a multieomponent droplet. The dot-dash eurve labeled “110% 
pure” is for a referenee ease with = 0 and Rf = 0.5 for pure butanol droplet at Pv,oo/p*sv,oo = 
1.10. The values of at Pv,oo/p*sv,oo = 0.6, 0.9, and 0.99 are 13, 172, and 8131. The values of at 
Pv,oo/p*sv,oo = 1-01 and 1.10 are 2488, and 514. For the dot-dash eurve with R^ = 0, the value of fg is 
649. 

values of are 2488 and 514, for oo/p*si;,oo = 1-01 and 1.10, respeetively. As a referenee, a dot-dash 
eurve is also shown in Fig. 3 for apure butanol droplet (with R^ = 0 and R^ = 0.5) !iipv,oo/p*sv,oo = 1.10, 
indieating slower growth of R{t) than a butanol-EG droplet espeeially near t* = 0. The value of tg for the 
dot-dash eurve with = 0 in Eig. 3 is 649. But with time, the dotted eurve will eonverge to the dot-dash 
eurve, beeause as the droplet grows indefinitely its volatile solvent mole fraetion will inerease toward unity 
asymptotieally. 

Similar to Eig. 3, the eondensation growth behavior of water-EG droplet is shown in Eig. 4. Beeause it 
is neeessary to have Pv,oo/p*sv,oo > 0.795 (= x(0), whieh is larger than that for the butanol-EG mixture) 
for eondensation growth to oeeur, a water-EG droplet is expeeted to grow at a slower rate than a butanol-EG 
droplet as refleeted in a eomparison of Eig.3 and Eig. 4. Eor the water-EG droplet with 7?^ = 0.4 and 
Rf = 0.5, the value of tggg is 2660 for Pv,oo/p*sv,oo = 0.99, and the values of are 4471 and 600, for 
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Figure 4: As in Fig. 3 but for water-EG droplet at Pv,oo/p*sv,oo = 0.99, 1.01, and 1.10. The value of t^gg at 
Pv,oo/p*sv,oo = 0.99 is 2660. The values of a.tpv^oc/p*sv,oo = 1-01 and 1.10 are 4471, and 600. For the 
dot-dash reference curve with Rn = 0, the value of is 649. 


Pv,oo/p*sv,oo = 1-01 and 1.10, respectively. The value of fg for the dot-dash curve with Rn = 0 in Fig. 
4 is 649 again (as independent of solvent type in the absence of nonvolatile cosolvent due to the definition 
of normalized time). A less significant difference between curves for the water-EG droplet (dotted) and 
pure water droplet (dot-dash) also suggests a weaker nonvolatile cosolvent effect on vapor transport for the 
water-EG droplet than for butanol-EG droplet, as noticed in the case of evaporating droplet. 

In general, condensation growth of droplet exhibits shallower slope than evaporating droplet, as illus¬ 
trated by comparing time scales shown in Eig. 1, Eig. 2 with that in Eig. 3, Eig.4. This may be explained by 
the fact that dR/dt ~ 1/i? according to ( [I^ , indicating a tendency of accelerating droplet size reduction 
during evaporation while decelerating droplet size growth during condensation. In the case of condensation 
growth, the volatile solvent mole fraction increases with time and therefore the nonvolatile cosolvent effect 
decreases progressively in the process. Thus, no physical mechanism exists for the acceleration-deceleration 
reversal phenomenon in condensational growth of a multicomponent droplet. While the presence of non¬ 
volatile cosolvents tends to hinder droplet evaporation, it consistently enhances the condensational growth 
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Figure 5: Surface temperature Tr evolution for water-EG droplet with nonvolatile cosolvent volume fraction 
= 0.4 and remnant volume fraction = 0.5, atp„^oo/p*sD,oo = 0 (solid), 0.2 (dot-dash), 0.99 (dashed), 
and 1.10 (dotted) according to (|2^. 


rate of a multicomponent droplet due to the effective reduction of saturation vapor pressure of the volatile 
solvent. 


3.3 Droplet Surface Temperature 


In theoretical analysis, an explicit expression of the local temperature at droplet surface is sometimes useful 
(e.g., Kulmala et al.[ 1993| Vesala et^ 19971. As a by-product of the mathematical derivation, the formula 
for local temperature at the multicomponent droplet surface Tr, according to Q and Q, can be obtained as 


Tr 


, npiLD dR^\ 

“ V ^ 2kT^ dt ) 


= T„ 


PlTT) 1 Pv,oo /sv,oo 
kT^ A + B/{R^-Rl)\ ’ 


(23) 


with dR^/dt, Tg, A, and B given by ( p^ — ( [T^ . Because B = 0 when = 0, droplet surface temperature 
Tr remains constant because dR?/dt = constant in the absence of the nonvolatile cosolvent. 
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For a water droplet {Rn = 0) at Pv,oo/p*sv,oo = 0, 0.2, 0.99, and 1.10 with Too = 300 K, its surface 
temperature Tr = 285.87, 288.69, 299.86, and 301.41 K, respectively. The presence of nonvolatile co¬ 
solvent leads to at time-dependent Tr for dR^/dt is not a constant anymore. Fig. 5 shows the curves of 
for water-EG droplet with R^ = 0.4 and R^ = 0.5 at Pv,oo/p*sv,oo = 0 (solid), 0.2 (dot-dash), 
0.99 (dashed), and 1.10 (dotted). For p„^oo/f*sd,oo < 1> each curve approaches the asymptotic value 300 K 
as dR^/dt —)• 0 with increasing t*, because the diminishing rate for vapor transport leads to an isothermal 
process. For Pv,oc/p*sv,oc = 1-10 (> 1) with R^ = 0.4 and Rf = 0.5 (dotted curve in Fig. 5), the value 
of Tr asymptotically approaches 301.41 K (> 300 K) due to a non-diminishing dR?/dt even as 72 —)■ cx) 
where dR?/dt —)• 2{pv^oo/p *sv,oo —1)This phemomanon is generally true for Pv,oo/p*sv,oo > 1 ac¬ 
cording to ( [T^ . The variation of surface temperature Tr appears to be much more dramatic for evaporating 
droplet than growing droplet by condensation. The “S-shape” curves in Fig. 5 for evaporating multicompo¬ 
nent droplet are consistent with the acceleration-deceleration reversal behavior exhibited in R{t*) curves in 
Fig. 1 and Fig. 2. Without the acceleration-deceleration reversal, the surface temperature evolution curves 
for a droplet during condensational growth appear less eventful. 

4 Summary 

Analytical formulas are derived in the present work for evaluating vapor transport of a volatile solvent for 
an isolated multicomponent droplet in a quiescent environment, based on quasi-steady-state approximation. 
Among multiple solvent components, only one component is considered to be much more volatile than 
the rest such that other components are assumed to remain unchanged in the droplet during the process of 
volatile solvent evaporation or condensation. With direct application to the Aerosol Jet® printing in mind, 
simplifications were made justifiably for ink droplef diamefer around a range of 1 fo 5 pm such as ignoring 
Ihe Knudsen number effecl, negligible vapor mass concenfrafion compared wifh fhe ambienf gas density, as 
well as applicability of Raoulf’s law for an ideal solufion of mulficomponenf solvenf mixfure. The form of 
derived equafion for dR?/dt suggesfs fhaf in fhe presence of nonvolatile cosolvenf, fhe rale of droplef size 
change diminishes as fhe amounl of volatile solvenf in fhe droplef approaches equilibrium value wifh fhe 
ambienf vapor pressure (72 —)• Re). In fhe absence of nonvolafile cosolvenf, fhe droplef radius 72 changes 
according lo fhe familier law”, which is recovered from fhe general formula derived herewifh when fhe 
volume fraclion of nonvolafile cosolvenf becomes zero (72„ = 0). 

In fhe case of evaporating droplef, fhe radius 72(7*) initially follows fhe “d? law” near 72(0) = 1 wifh 
an accelerafed rale of change. However, fhe presence of nonvolafile cosolvenf is predicled fo conlinuously 
reduce fhe rale of evaporafion as fhe mole fraclion of fhe volafile solvenf keeps decreasing. According fo 
Raoull’s law, decreasing fhe mole fraclion effeclively corresponds lo a reducfion of salurafion vapor pressure 
for fhe volafile solvenf and fhereby weakens fhe driving force for evaporafion. Thus, fhe magnilude of droplef 
size change is evenlually decelerated toward fhe end wifh dR/dt —)■ 0 as 72 —)■ Re- Such an acceleralion- 
deceleralion reversal behavior in fhe droplef size change is unique in fhe presence of nonvolatile cosolvenf, 
white fhe droplef of pure solvenf follows fhe d? law of accelerating size reducfion all fhe way Ihrough fhe 
end. A closer examinalion of fhe explicif formula for d'^R/dt^, however, indicales fhaf fhe acceleralion 
phase may disappear when fhe nonvolafile cosolvenf volume fraclion 72^ is small and fhe ambienf vapor 
pressure Pv,oo is relafively high. 

Because fhe nef effecl of adding nonvolafile cosolvenf is to reduce fhe mote fraction of fhe volatile 
solvenf such fhaf fhe salurafion vapor pressure is lowered, vapor condensalion onto fhe mulficomponenf 
droplef is predicled fo occur when fhe ambienf vapor pressure is subsalurafed wifh respecf lo fhaf for fhe 
pure volafile solvenf. In Ihis case, fhe droplef will grow asympfolically loward a finite size 72 = 72e (> 1). 
Such resfricled growlh of droplef by condensalion is unique only when fhe nonvolatile cosolvenf is presenl 
in fhe mulficomponenf droplef. Buf when fhe ambienf vapor pressure becomes supersaluraled wifh respecf 
fo fhaf for fhe pure volafile solvenf, fhe condensation growlh of droplef can conlinue indefinitely wifhoul 
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bound. Condensation growth of droplet generally exhibits shallower slope than evaporating droplet beeause 
the droplet size ehange is eontinuously deeelerated with time as illustrated in figures of R{t*) and 

The effeets of nonvolatile eosolvent on the volatile solvent evaporation-eondensation eharaeteristies 
predieted here should oeeur with general multieomponent droplets, in view of the faet that most multieom- 
ponent mixtures eontain eosolvents with different volatilities. It is expeeted more often than not to observe 
a single eomponent of the most volatile solvent evaporating or eondensing at a mueh shorter time seale than 
the other eosolvents. Thus, the formulas derived here may beeome applieable to more general situations 
with multieomponent droplets as long as the evaporation-eondensation time seales for different eomponents 
differ substantially. 


Acknowledgments 

The author would like to thank John Lees for eneouragement and guidanee, and Dr. Mike Renn, Dr. Kurt 
Christenson, Jason Paulsen, as well as many other Optomee eolleagues, for helpful teehnieal diseussions. 


References 

Aggarwal, S. K. (1998). A review of spray ignition phenomena: Present status and future researeh. Prog. 
Energy Combust. Sei. 24:565-600 

Bird, R. B., Stewart, W. E., and Eightfoot, E. N. (1960). Transport Phyenomena. John Wiley, New York 

Christenson, K. K., Paulsen, J.A., Renn, M.J., MeDonald, K., and Bourassa, J. (2011). Direet printing of 
eireuit boards using Aerosol Jet®. Proe. NIP 27 Digital Eabrie. 433^36 

Davies, J. E, Haddrell, A. E., and Reid, J. P. (2012). Time-resolved measurements of the evaporation of 
volatile eomponents from single aerosol droplets. Aerosol Sei. Teehnol. 46(6):666-677 

Peng, J. Q. (2013). Diffusion-eontrolled quasi-stationary mass transfer for an isolated spherieal partiele in 
an unbound medium. Chem. Eng. Comm. 200:1:65-76 

Kahn, B. E. (2007). The M^D aerosol jet system, an alternative to inkjet printing for printed eleetronies. 
Organie and Printed Eleetronies. 1:14-17 

Kulmala, M., Vesala, T., and Wagner, PE. (1993). An analytieal expression for the rate of binary eondensa- 
tional partiele growth. Proe. R. Soe. Eond. A 441:589-601 

Pin, J.-C. and Gentry, J. W. (2003). Spray drying drop morphology: Experimental study. Aerosol Sei. Teeh¬ 
nol. 37:15-32 

Maxwell, J. C. (1878). Diffusion, in Eneyelopaedia Britanniea, Vol. 7 (9th ed.), Eondon, 214-221 

Peng, C., Chow, A. H. E., and Chan, C. K. (2000). Study of the hygroseopie properties of seleeted pharma- 
eeutieal aerosols using single partiele levitation. Pharm. Res. 17:1104-1109 

Pruppaeher, H. R. and Klett, J. D. (1978). Mierophysies of Clouds and Preeipitation. D. Reidel Publishing, 
Dordreeht 

Ravindran, P. and Davis, E. J. (1982). Multieomponent evaporation of single aerosol droplets. J. Colloid 
Interf. Sei. 85:278-288 

Seriven, E. E. (1959). On the dynamies of phase growth. Chem. Eng. Sei. 10:1-13 


15 



Vesala, T., Kulmala, M., Rudolf, R., Vrtala, A., and Wagner, R E. (1997). Models for condensational growth 
and evaporation of binary aerosol particles.. J. Aerosol Sci. 28(4):565-598 

Widmann, J. F. and Davis, E. J. (1997). Evaporation of multicomponent droplets. Aerosol Sci. Technol. 
27:243-254 

Wilms, J. (2005). Evaporation of multicomponent droplets. PhD thesis. University of Stuttgart 

Zollmer, V., Muller, M., Renn, M., Busse, M., Wirth, I., Codlinski, D., Kardos, M. (2006). Printing with 
aerosols: A maskless deposition technique allows high definition printing of a variety of functional mate¬ 
rials. Euro. Coating J. 07-08:46-55 


Figure Captions 


Figure 1: Evaporating droplet radius evolution with time for water-EG and butanol-EG mixtures, with 
ethylene glycol (EG) serving as the nonvolatile cosolvent with a volume fraction of = 0.4. The remnant 
volume of droplet (excluding that of the volatile solvent) is denoted as 47r7?^/3 here in this figure with the 
remnant volume fraction Rf = 0.5. The solid curve and dashed curve correspond to Pv,oo/p*sv,oo = 0 and 
0.2, with = 1.23 and 1.71 for water-EG droplet while = 4.81 and 5.49 for butanol-EG droplet. 
The dotted curve is for pure volatile solvent with = 0 and Rf. = 0.5 at Pv,oo/p*sv,oo = 0, as a reference. 

Figure 2: As in Fig. 1 but for R^ = 0.1 and Rf = 0.2, at Pv,oo/p*sv,oo = 0 and 0.2, with = 1.06 and 
1.37 for water-EG droplet while = 1.94 and 2.57 for butanol-EG droplet. 

Figure 3: As in Fig. 1 (with nonvolatile cosolvent volume fraction = 0.4 and remnant volume fraction 
Rf = 0.5) but only for butanol-EG droplet at Pv,oo/p*sv,oo = 0.6, 0.9, 0.99, 1.01, and 1.10 to illustrate 
characteristics of condensation growth of a multicomponent droplet. The dot-dash curve labeled “110% 
pure” is for a reference case with R^ = 0 and Rf = 0.5 for pure butanol droplet at Pv,oo/p*sv,oo = 
1.10. The values of fggg at Pv,oo/p*sv,oo = 0.6, 0.9, and 0.99 are 13, 172, and 8131. The values of at 
Pv,oo/p*sv,oo = 1-01 and 1.10 are 2488, and 514. For the dot-dash curve with R^ = 0, the value of fg is 
649. 


Figure 4: As in Fig. 3 but for water-EG droplet at Pv,oo/p*sv,oo = 0.99, 1.01, and 1.10. The value of tggg at 
Pv,oo/p*sv,oo = 0.99 is 2660. The values of fg at Pv,oo/p*sv,oo = 1-01 and 1.10 are 4471, and 600. For the 
dot-dash reference curve with Rn = 0, the value of fg is 649. 


Figure 5: Surface temperature Tr evolution for water-EG droplet with nonvolatile cosolvent volume fraction 
R^ = 0.4 and remnant volume fraction = 0.5, atpt,,c!o/p*si;,oo = 0 (solid), 0.2 (dot-dash), 0.99 (dashed), 
and 1.10 (dotted) according to (231. 
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